Error estimation and adaptivity for incompressible, non–linear (hyper–)elasticity

A Galerkin finite element method is developed for non–linear, incompressible (hyper) elasticity, and a posteriori error estimates are derived for both linear functionals of the solution and linear functionals of the stress on a boundary where Dirichlet boundary conditions are applied. A second, higher order method for calculating a linear functional of the stress on a Dirichlet boundary is also presented together with an a posteriori error estimator for this approach. An implementation for a 2D model problem with known solution demonstrates the accuracy of the error estimators. Finally the a posteriori error estimate is shown to provide a basis for effective mesh adaptivity

A Comparison of Numerical Methods used for\ud Finite Element Modelling of Soft Tissue\ud Deformation

Soft tissue deformation is often modelled using incompressible nonlinear elasticity, with solutions computed using the finite element method. There are a range of options available when using the finite element method, in particular, the polynomial degree of the basis functions used for interpolating position and pressure, and the type of element making up the mesh. We investigate the effect of these choices on the accuracy of the computed solution, using a selection of model problems motivated by typical deformations seen in soft tissue modelling. We set up model problems with discontinuous material properties (as is the case for the breast), steeply changing gradients in the body force (as found in contracting cardiac tissue), and discontinuous first derivatives in the solution at the boundary, caused by a discontinuous applied force (as in the breast during mammography). We find that the choice of pressure basis functions are vital in the presence of a material interface, higher-order schemes do not perform as well as may be expected when there are sharp gradients, and in general that it is important to take the expected regularity of the solution into account when choosing a numerical scheme

Cardiac Electromechanics: The effect of contraction model on the mathematical problem and accuracy of the numerical scheme

Models of cardiac electromechanics usually contain a contraction model determining the active tension induced at the cellular level, and the equations of nonlinear elasticity to determine tissue deformation in response to this active tension. All contraction models are dependent on cardiac electro-physiology, but can also be dependent on\ud the stretch and stretch-rate in the fibre direction. This fundamentally affects the mathematical problem being solved, through classification of the governing PDEs, which affects numerical schemes that can be used to solve the governing equations. We categorise contraction models into three types, and for each consider questions such as classification and the most appropriate choice from two numerical methods (the explicit and implicit schemes). In terms of mathematical classification, we consider the question of strong ellipticity of the total strain energy (important for precluding ‘unnatural’ material behaviour) for stretch-rate-independent contraction models; whereas for stretch-rate-dependent contraction models we introduce a corresponding third-order problem and explain how certain choices of boundary condition could lead to constraints on allowable initial condition. In terms of suitable numerical methods, we show that an explicit approach (where the contraction model is integrated in the timestep prior to the bulk deformation being computed) is: (i) appropriate for stretch-independent contraction models; (ii) only conditionally-stable, with the stability criterion independent of timestep, for contractions models which just depend on stretch (but not stretch-rate), and (iii) inappropriate for stretch-rate-dependent models

Party membership and campaign activity in Britain: The impact of electoral performance

The article examines the impact of electoral results on party membership and activity. Previous studies have focused on the long-term effects of electoral success or failure, suggesting that they may produce a spiral of demobilization or mobilization. The article shows that the dramatic change of electoral fortunes experienced by British parties at the 1997 general election broke this spiral, with the outcome leading to significant changes in the health and activity of local parties. It is concluded that dramatic election results can have significant implications for party organization

Error bounds on block Gauss Seidel solutions of coupled\ud multiphysics problems

Mathematical models in many fields often consist of coupled sub–models, each of which describe a different physical process. For many applications, the quantity of interest from these models may be written as a linear functional of the solution to the governing equations. Mature numerical solution techniques for the individual sub–models often exist. Rather than derive a numerical solution technique for the full coupled model, it is therefore natural to investigate whether these techniques may be used by coupling in a block Gauss–Seidel fashion. In this study, we derive two a posteriori bounds for such linear functionals. These bounds may be used on each Gauss–Seidel iteration to estimate the error in the linear functional computed using the single physics solvers, without actually solving the full, coupled problem. We demonstrate the use of the bound first by using a model problem from linear algebra, and then a linear ordinary differential equation example. We then investigate the effectiveness of the bound using a non–linear coupled fluid–temperature problem. One of the bounds derived is very sharp for most linear functionals considered, allowing us to predict very accurately when to terminate our block Gauss–Seidel iteration.\ud \ud Copyright c 2000 John Wiley & Sons, Ltd

Nonuniqueness in a minimal model for cell motility

Two–phase flow models have been used previously to model cell motility, however these have rapidly become very complicated, including many physical processes, and are opaque. Here we demonstrate that even the simplest one–dimensional, two–phase, poroviscous, reactive flow model displays a number of behaviours relevant to cell crawling. We present stability analyses that show that an asymmetric perturbation is required to cause a spatially uniform, stationary strip of cytoplasm to move, which is relevant to cell polarization. Our numerical simulations identify qualitatively distinct families of travelling–wave solution that co–exist at certain parameter values. Within each family, the crawling speed of the strip has a bell–shaped dependence on the adhesion strength. The model captures the experimentally observed behaviour that cells crawl quickest at intermediate adhesion strengths, when the substrate is neither too sticky nor too slippy

A note on heat and mass transfer from a sphere in Stokes\ud flow at low Péclet number

We consider the low Péclet number, Pe ≪ 1, asymptotic solution for steady-state heat and mass transfer from a sphere immersed in Stokes flow with a Robin boundary condition on its surface, representing Newton cooling or a first-order chemical reaction. The application of van Dyke’s rule up to terms of O(Pe3) shows that the O(Pe3 log Pe) terms in the expression for the average Nusselt/Sherwood number are double those previously derived in the literature. Inclusion of the O(Pe3) terms is shown to increase significantly the range of validity of the expansion

Political participation: the vocational motivations of Labour party employees

Party employees are an under-researched group in political science. This article begins to address this oversight by examining Labour Party employees using new quantitative and qualitative data. It argues that party employment should be regarded as a form of political participation and as a consequence, existing models of political participation can be utilised to help explain why people work for political parties. After testing these propositions, the article concludes that existing models are indeed helpful in explaining the motivations for party employment

On the deflection of a liquid jet by an air-cushioning layer

A hierarchy of models is formulated for the deflection of a thin two-dimensional liquid jet as it passes over a thin air-cushioning layer above a rigid flat impermeable substrate. We perform a systematic derivation of the leading-order equations of motion for the jet in the distinguished limit in which the air pressure jump, surface tension and gravity affect the displacement of the centreline of the jet, but not its thickness or velocity. We identify thereby the axial length scales for centreline deflection in regimes in which the air layer is dominated by viscous or inertial effects. The derived length scales and reduced equations aim to expand the suite of tools available for future analyses of the evolution of lamellae and ejecta in impact problems. Assuming that the jet is sufficiently long that tip and entry effects can be neglected, we demonstrate that the centreline of a constant-thickness jet moving with constant axial speed is destabilised by the air layer for sufficiently small surface tension. Expressions for the fastest-growing modes are obtained in both the viscous-dominated air and inertia-dominated air regimes. For a finite-length jet emanating from a nozzle, we show that, in one particular asymptotic limit, the evolution of the jet centreline is akin to the flapping of an unfurling flag above a thin air layer. We discuss the distinguished limit in which tip retraction can be neglected and perform numerical investigations into the resulting model. We show that the cushioning layer causes the jet centreline to bend, leading to rupture of the air layer. We discuss how our toolbox of models can be adapted and utilised in the context of recent experimental and numerical studies of splash dynamics